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Abstract 

The superconducting proximity effect is investigated for SN double 
layers in a regime where the resulting transition temperature T c does 
not depend on the mean free paths of the films and, within limits, not 
on the transparency of the interface. This regime includes the thin 
film limit and the normalized initial slope S sn = {d s /T s )\dT c / dd n \. 
The experimental results for T c are compared with a numerical sim- 
ulation which was recently developed in our group. The results for 
the SN double layers can be devided into three groups: (i) When N 
= Cu, Ag, Au, Mg a disagreement between experiment and theory by 
a factor of the order of three is observed, (ii) When N = Cd, Zn, Al 
the disagreement between experiment and theory is reduced to a fac- 
tor of about 1.5, (iii) When N = In, Sn a reasonably good agreement 
between experiment and theory is observed. 
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1 Introduction 



The properties of a superconducting film or thin wire S are modified when 
they are in contact with a normal metal N. This phenomenon was first ob- 
served in the pioneering experiments by Meissner who explored the prop- 
erties of superconducting wires covered with normal metals. It is generally 
called the "superconducting proximity effect" (SPE). It was intensively stud- 
ied in the 1960 's [2], [3], [I], [3], 0, During the last decade it has expe- 
rienced a renewed interest theoretically g], 0, HQ], HH, [12], P3], d, 113, 
DH, D3, GB], d, EH as well as experimentally EH, E2, E3, El, [IS], 
PH, E3, EH], E3, EE]- Recently the SPE has been extended to SN-multi 
layers EH], [HI- 

A few years ago, our group EZj investigated the proximity effect between 
Pb and several alkali metals. For a better analysis of these measurements 
we developed a quantitative numerical method for the calculation of the 
transition temperature of an SN double layer [31]. Our numerical results 
show that, when a superconductor S is covered with a normal metal N, that 
the initial slope ^p- is independent of the mean free paths of the two metals 
and the transparency of the interface (if the transmission is not dramatically 



changed). If one defines a normalized initial slope S sn — ' 



then S* 



T 3 

is independent of the thickness d s of the superconductor up to relatively 
large values of d s . If the superconductor is very weak coupling (2irT s « 
©d, Oz?=Debye temperature) then our result for the initial slope converges 
towards the results for the thin film Cooper limit [32] (see below). 

When we compared our experimental initial slope ^| dn=0 with our nu- 
merical calculation we observed that the experimental results were consider- 
ably smaller than the theoretical predictions. Surprised by the discrepancy 
we searched the literature and found early experiments from the 1960s, par- 
ticularly by Hilsch [2], [H] and Minigerode [7], from which the normalized 
initial slope can be derived. These measurements showed a similar disagree- 
ment in the initial slope with the theory (see ref. p?T]). 

Since we were rather amazed by the discrepancy between our experiments 
and theory in the SPE and also by the fact that this discrepancy had not 
been detected previously, we decided to re-investigate the SPE. In this paper 
we investigate the SPE in the range, where a minimum of experimental pa- 
rameters is needed to perform a quantitative comparison with the theory. We 
focus on the normalized initial slope S sn of SN sandwiches and the transition 
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temperature of very thin NS sandwiches in the thin film limit. 



2 Experiment and Results 



We use thermal evaporation to condense the thin films onto a substrate at 
liquid helium temperature. To obtain clean films all the evaporation sources 
are surrounded with liquid N 2 and the vacuum in our system is better than 
10 _11 torr. We first evaporate 10 atomic layers of insulating Sb on a helium 
cold quartz substrate. The Sb film acts as a fresh substrate and insures that 
the following quench condensed films are flat and homogeneous. 

In a series of experiments a film of the superconductor Pb is first con- 
densed onto the Sb substrate. Afterwards the Pb is covered in several step 
with an increasing thickness of the normal metal Ag. The thickness of the 
films is measured with a quartz oscillator. The accuracy of the thickness 
measurement is about 15%. After each evaporation the superconducting 
transition curve, the magnetoresistance and the Hall effect of the double 
layer are measured. Fig.l shows a plot of T c versus the Ag thickness d Ag on 
top of a 251A thick Pb film. This plot yields graphically the initial slope 
dT c /dd n \d n =o and the normalized initial slope S sn 



where dpb = d s is the thickness of the Pb films and T s = 7.2K is the transition 
temperature for Pb alone. 




4 dT c (d n = 0) 
T s dd n 



(1) 
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Fig.l: T c versus d\ g for an PbAg double layer. 

This experiment is repeated for different thicknesses of the superconduc- 
tor Pb. In Fig.2 S sn is plotted versus the Pb thickness dpb- It is essentially 
independent of the Pb thickness. This was the prediction of our numerical 
results. The value of the normalized initial slope is Spf,A g = 0.66 ± 0.05. 
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Fig.2: The normalized initial slope S sn for the Pb/Ag 
sandwiches as a function of the Pb thickness 



In the next experimental series we investigate NS double layers in the thin 
film limit and express the results in terms of the normalized initial slope S sn . 
As an example a thin Ag film (d,A g = 41. OA) is condensed onto the insulating 
Sb substrate and covered in several steps with increasing thickness of Pb. 
Fig. 3a shows the results. The inverse T c -reduction 1/AT C = 1/ (T s — T c ) is 
plotted as a function of the Pb thickness. From these plots we extract the 
value of the normalized initial slope. For the AgPb double layer the value is 
SpbAg = 0.66. This is in excellent agreement with the results from the first 
experimental series. 
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Figure 3. The inverse T c reduction 1/ (T s — T c ) versus the Pb thickness dpb 
for double layers of N/Pb where N stands for the metals Ag, Mg, Zn and Sn. 



We use this much more efficient procedure to investigate sandwiches of 
Pb with the normal metals Ag, Au, Cu, Mg. Furthermore we included also as 
"normal metals" superconductors with a transition temperatures T n which 
lie below the value of T c for Pb. These metals are Zn, Cd, Al, In and Sn. 
Fig.3b-d shows some of the results for Mg, Zn and Sn. 
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In table I the experimental data are collected. The first column gives 
the experimental code; which includes the symbols for the normal metal and 
the superconductor. The following columns give the thickness of the normal 
metal (or the superconductor with the lower T c ), the transition temperature 
T n of N if N is a superconductor, the ratio of the experimental (i.e., phonon 
enhanced) density of states and the experimental normalized initial slope 
S exp . The estimated error of S cxp is about 15%. The density of states are 
taken from Kittel's book [33]. The theoretical prediction of the normalized 
initial slope by the Werthamer theory Swh and our numerical result S S i m will 
be discussed below. The transition temperatures for quench condensed Cd 
and Zn are taken from ref. and j^Hj, the transition temperatures of 

quench condensed Al, In and Sn are taken from and [3*7j . 



exp. 


d n (A) 


T n (K) 


n;/n; 


c 

'-'exp 


Swh. 


c . 

kJ sim. 


^cxp/^sim. 


W (n;/n:) 


AgPbJB 


41.2 





0.387 


0.66 


0.95 


1.74 


0.38 


4.5 


CuPbJE 


32.9 





0.603 


0.98 


1.49 


2.67 


0.37 


4.43 


AuPbJD 


29.5 





0.442 


0.62 


1.09 


1.98 


0.31 


4.48 


MgPbJG 


28.5 





0.572 


0.48 


1.41 


2.54 


0.19 


4.44 


CdPbJJ 


31.4 


0.80 


0.329 


0.39 


0.40 


0.528 


0.74 


1.60 


ZnPbUU 


26.0 


1.39 


0.430 


0.33 


0.44 


0.54 


0.61 


1.26 


AlPbUT 


21.5 


2.28 


0.833 


0.44 


0.68 


0.76 


0.58 


0.91 


InPbUP 


32.1 


4.1 


0.663 


0.30 


0.31 


0.354 


0.85 


0.53 


SnPbUN 


33.0 


4.7 


0.664 


0.29 


0.25 


0.283 


1.02 


0.43 





















Table I: The experimental normalized initial slope S exp are compared with 
the theoretical predictions by Werthamer's theory 5*wh an d the author's 
numerical calculations S S i m . The first four columns give the experimental 
code (containing the symbols of the normal conductor N and superconductor S), 
the thickness of N, the transition temperature of N (if superconducting) and 
the ratio of the experimental density of states N*/N* (which includes 
electron-phonon enhancement). The second to last column gives the ratio 
between the experimental S cxp and the numerical results S S i m - 

3 Discussion 

Let us consider a double layer with the superconducting transition tempera- 
ture T c . This transition temperature defines a characteristic time tt c = t c = 
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h/ (irk B T c ). The superconducting coherence length is the distance that an 
electron propagates during this time t c . In a clean metal this is equal to 
£o — vft c and in a dirty metal one has £ = \/Dt c where D is the diffusion 
constant of the dirty metal. 

When, during the time r c , an arbitrary conduction electron propagates 
through the whole thickness range of the double layer with roughly equal 
probability then the superconducting properties of the system are averaged 
over both films. In this case the mean free paths of the individual metal have 
no bearing on the superconducting transition temperature T c . (However, the 
mean free paths determine whether the system is in this limit). A well-known 
example is the thin film or " Cooper" limit of a double layer when both films 
are much thinner than their coherence lengths and the "transmission time" 
through the interface is much shorter than r c In the limit that T c is much 
smaller than the Debye temperature, T c << O^, the transition temperature 
of such a double layer is given by the BCS-Cooper formula 



T c = 1.146b exp 



(NV) ef/ 

where the effective BCS interaction (NV) e j of the double layer is 

d s N s (NV) a + d n N n (NV) n 



(2) 



(NV)ef 



d s N s + d n N n 



and (NV) sn are the BCS interactions in the super- and normal conductor. 

Using the BCS-Cooper formula for an SN double layer yields for the 
normalized initial slope 



c d s 

s CP - ¥ 



dT n 



ddn 



d. 



d(\n(T c )) 



ddn 



N n (NV) S - (NV), 



((NV) s y 



(3) 



The predictions of equation (jHJ) for the normalized initial slope (which is 
derived from equation ©) has two problems, (i) Pb is not a superconductor 
with T c << Q D and therefore the equation © does not yield a good repre- 
sentation of Pb, (ii) equation (J2J) uses the same Debye temperature for both 
metals and this is generally not fulfilled in the experiment. One can replace 
the BCS-Cooper formula for T c by an (implicit) expression 



E 



(4) 
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where n c = Qp/ (2tiT c ). We do not evaluate equation (JIJ) here because its 
results are in included in our numerical procedure as a limiting case. 

Werthamer jHS] derived a set of equations for the transition temperature 
of a double layer of a super- and a normal conductor in the dirty limit, using 
de Gennes' interface boundary condition In this theory the gap pa- 

rameter is approximated as a cosine and hyperbolic cosine function in the 
superconductor and normal conductor respectively. (This is sometimes called 
the single mode expansion). This theory agrees quite well with the experi- 
mental results for double layers of two superconductors The normalized 
initial slope can be expressed as 

5wh = ^T^ 1( - ln( | )} (5) 

here x -1 (y) is the inverse function of x{ x ) = + \ x ) ~ ^(f) an d ^( z ) is 
the digamma function. If we assume that the transition temperature T n for 
the normal metals Cu, Ag, Au and Mg is infinitely small, then x~ l {~ ^(t 8 ")) 
takes the value one. The values of Swh according to Werthamer's theory are 
included in table I. The values of S"wh don't show a good agreement for the 
normal metals Cu, Ag, Au and Mg. 

Finally we compare the experimental values 5 exp with our numerical 
calculation. This numerical calculation derives the transition temperature 
of a double or multi layer of a superconductor and a normal conductor. 
The equivalence in the propagation of the superconducting pair amplitude 
and a single electron in Gorkov's linear gap equation is used. The single 
electrons act as messengers which carry the information about the super- 
conducting gap (N S A (r') /tt) from one position-time (r', t' < 0) to another 
position-time (r, t = 0). This message which decays thermally with time as 
Vt (t) = Y^\u) n \<n D ex P ( — ^ \ Un \ \t'\)-> ls integrated at (r,t = 0) over all starting 
position-times (r' ; t') and, after multiplication with the BCS interaction V s , 
yields the new gap function A (r). At the transition temperature the proce- 
dure has to be self-consistent, i.e. the initial and final gap function have to 
be identical. The propagation of the single electrons is then quasi-classically 
simulated. The frame work of the calculation is the weak coupling theory of 
superconductivity. 

This numerical procedure to calculate the transition temperature of dou- 
ble or multi-layers consisting of thin films of superconductors and normal 
conductors is very flexible. It uses the following input parameters of the 
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individual metal films (i) thickness, (ii) density of states and Fermi velocity, 
(iii) transition temperature, (iv) Debye temperature, (v) mean free path and 
(vi) transmission through the interface between the films. 

An important outcome of the numerical simulation is the result that the 
normalized initial slope of an SN double layer as a function of d n at d n = 
does not depend on (i) the mean free paths of the two metals, (ii) the 
thickness of the superconductor and (iii) a finite (but not too large) barrier 
between the two metals. 

We include the numerical results S S i m in table I. For the double layers with 
the normal conducting metals Cu, Ag, Au and Mg the discrepancy between 
experiment and theory is very large, of the order of 3. It is remarkable that 
the deviation is considerably smaller when the "normal metal" is really a 
superconductor with a smaller transition temperature although the deviation 
is still factor of about 1.5. For In and Sn, however, there is a good agreement 
between theory and experiment. 

It is surprising that this disagreement between experiment and theory 
has not been noticed in the past. The main reason is that the majority of 
the experimental and the theoretical work focused on NS sandwiches with 
thick normal metal films. Then superconductivity is only obtained for a 
finite thickness of the superconductor. In this case a comparison between 
experiment and theory requires many fit parameters such as the transparency 
of the interface and the mean free paths of the superconductor and the normal 
conductor. Therefore it is quite possible to fit the experimental data by using 
the wrong parameters that can't be checked otherwise. 

The physical origin of this disagreement between experiment and theory 
is not understood. Our theoretical simulation of the SPE uses the frame 
work of weak coupling superconductivity. Quench condensed Pb, In and Sn 
are not weak coupling. The ratios of 2A / (ksTg) for quench condensed films 
are 4.6 for Pb, 3.9 for In and 4.0 for Sn j^Zj, jlU]. An obvious proposal 
would be to solve the superconducting proximity effect for strong coupling 
superconductors. This means to develop and solve a series of equations for 
the energy and position dependent gap function A (r, uj{). This would be 
an extremely demanding job. As a start, we considered the thin film limit 
for the SPE of strong superconductors. This consideration, which will be 
published elsewhere, does not remove the discrepancy between experiment 
and theory. 

The fact that the Werthamer result 5*wh disagrees less with the experi- 
mental data S exp for Cu, Ag, Au and Mg than our numerical result might be 
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accidental. The single mode expansion of A (r) in the Werthamer theory is 
surely less appropriate than a self-consistant gap function. 

The experimental normalized initial slope is proportional to the density 
of states ratio. Although the density of states can be modified in quench 
condensed films it is inconceivable that this explains a factor of three in the 
initial slope. Let us return to Fig.l where the reduction of T c of Pb by a thin 
layer of Ag is plotted. 

4 Conclusion 

In this paper, the superconducting proximity effect is investigated for SN 
double layers in a regime where the resulting T c does not depend on the 
mean free path of the films and, within limits, not on the transparency of 
the interface. This regime includes the thin film limit and the normalized 
initial slope S sn = (d s /T s )\dT c /dd n \ . While the layer S is always Pb the 
layer N is either a non-superconducting metal such as Cu, Ag, Au and Mg 
or a superconductor with a T c below the transition temperature of Pb. The 
experimental results for the transition temperature T c are compared with a 
numerical calculation which was recently developed in our group. The results 
for the SN double layers can be divided into three groups: 

• When N represents a non-superconducting metal film (N=Cu, Ag, Au 
and Mg) we observe grave deviations between experiment and theory 
by a factor of the order of three. 

• When N represents a superconductor with a low T c (N=Cd, Zn, Al) 
the deviation between experiment and theory is still there but reduced 
by a factor of two. 

• When N represents a superconductor with a T c which is about half 
the T c of Pb (N=In, Sn) then we observe a reasonably good agreement 
between experiment and theory. 

Prior to our recent experiments we believed that the proximity effect be- 
tween a superconductor and a normal conductor represented an intensively 
studied phenomenon with a good theoretical understanding. We are deeply 
puzzled by the large observed discrepancy between experiment and theory. 
It would be very desirable if other theoretical approaches would give quanti- 
tative predictions for the normalized initial slope in SN double layers. There 
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have been a number of theoretical papers published which extended the prox- 
imity effect to more complex systems, for example between a superconductor 
and a ferromagnet but which include implicitly the simpler case of an SN 
double layer. These authors would be able to calculate quantitatively the 
normalized initial slope from their theory. 

Experimentally it would be desirable to extend the measurements to SN 
layers where S is a weak coupling superconductor. This requires lower tem- 
peratures but permits the use of thicker films because the coherence lengths 
are larger at lower temperatures. Aluminum would be a good candidate for 
the superconductor if evaporated in ultra-high vacuum so that the Al is not 
granular. 
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